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Abstract. We simulate the dynamics of Au+Au collisions at the Relativistic 
Heavy Ion Collider (RHIC) with a hybrid model that treats the quark-gluon plasma 
macroscopically as an ideal fluid, but models the hadron resonance gas microscopically 
using a hadronic cascade. We find that much of the mass-ordering pattern for V2(pt) 
observed at RHIC is generated during the hadronic stage due to build-up of additional 
radial flow. We also find that the mass-ordering pattern is violated for <f> meson due 
to small interaction cross section in the hadron resonance gas. 



1. Introduction 

Whether the quark-gluon plasma (QGP) behaves like a "perfect liquid" is one of the 
important question in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) 
[H 121 [3J H]. The observed elliptic flow parameter V2 and its transverse momentum 
dependence agree well with predictions from ideal fluid dynamics assuming zero viscosity 
[51 E]- The ideal fluid dynamical description, however, gradually breaks down as one 
studies peripheral collisions or moves away from midrapidity [6] . This requires a more 
realistic treatment of the early and late stages in dynamical modeling of relativistic 
heavy ion collisions. In previous work [?1 [8] we have shown that a large fraction of 
these deviations from ideal hydrodynamics is due to "late viscosity". Here we report 
additional results from the hybrid model study focusing our attention on a detailed 
investigation of dissipative effects during the late hadronic rescattering stage. 

2. Model 

For the space-time evolution of the perfect QGP fluid we solve numerically the equations 
of motion of ideal fluid dynamics for a given initial state in three spatial dimensions 
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and in time [HJ [TU]. For the high temperature (T > T c = 170 MeV) QGP phase we 
use the equation of state of massless parton gas (u, d, s quarks and gluons) with a 
bag pressure B. We switch from ideal hydrodynamics to a hadronic cascade model 
at the switching temperature T sw = 169 MeV. The subsequent hadronic rescattering 
cascade is modeled by JAM [UJ, initialized with hadrons distributed according to the 
hydrodynamic model output, calculated with the Cooper- Frye formula [12] along the 
T sw = 169 MeV hypersurface. JAM implements experimental hadronic scattering cross 
section data where available and uses the additive quark model where data do not exist, 
assuming the following formula for the total cross section [TTJ [13[ [HI [15] : 

Here crj^y is the total nucleon-nucleon cross section, rij is the number of constituent 
quarks in a hadron and n S i is the number of strange quarks in a hadron. For hadrons 
composed entirely of strange quarks, such as (f> = (ss) and fl = (sss), the cross sections 
become very small due to the suppression factors in brackets in Eq. (TjQ). We note that, 
to study <fi mesons in our hybrid model, we stabilize them by turning off their decay 
channels during the hadronic cascade. For initial conditions in hydrodynamic equations, 
we employ the Glauber model suitably generalized to account for the longitudinal 
structure of particle multiplicity [7[ [16] . 



3. Results 




Figure 1. Transverse momentum dependence of the elliptic flow parameter «2(j»t)- 
(Left) i^Cpt) for pions, kaons and protons from the hybrid model are compared with 
STAR data [T7]. (Right) Solid (dashed) lines are with (without) hadronic rescattering. 



In Fig. 1 (left), we compare i>2(j>t) for pions, kaons and protons from the hybrid 
model with the STAR data [17J. We reasonably reproduce mass-splitting behaviour seen 
in the data. We note that we also reproduce the data in other centralities (not shown) 
except for very central collisions due to absence of eccentricity fluctuation [18]. Figure 
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1 (right) shows how mass-splitting of p^pr) is generated during evolution by switching 
off hadronic rescatterings in a hadronic cascade. Slope of pion ^(pr) becomes steeper 
due to additional development of elliptic flow and reduction of mean p^ |19j . For heavy 
hadrons, on the other hand, radial flow reduces at low pp [20]. Assuming positive 
elliptic flow, v±(tp—0, ir) > v± (</?=§ > T")> the stronger transverse flow v± in the reaction 
plane pushes heavy particles to larger p^ more efficiently in the reaction plane than 
perpendicular to it. The generation of additional radial flow in the hadronic stage is 
responsible for the mass-splitting of i>2(pr) observed in the low px region. From these 
observations we conclude that the large magnitude of the integrated vi and the strong 
mass ordering of ^(pr) observed at RHIC result from a subtle interplay between perfect 
fluid dynamics of the early QGP stage and dissipative dynamics of the late hadronic 
stage. 




Figure 2. (Left) U2(pt) from the hybrid model. (Right) «2(pt) from ideal 
hydrodynamics with T dec =100 MeV. Solid, dashed and dotted lines are results for 
</> mesons, pions and protons. 



(J) mesons have considerably smaller scattering cross sections than non-strange 
hadrons [21]. They are therefore expected to show larger dissipative effects in our 
hybrid model and not to fully participate in the additional radial flow generated during 
the hadronic rescattering stage. Figure 2 shows ^(pr) from the hybrid model (left) 
and the ideal hydrodynamics (right) for tt, p and 4>. As a result of rescattering the 
proton elliptic flow ends up being smaller than that of the (f> meson, vf(pr) < v t{PT) for 
<pr < 1-2 GeV/c, even though m ( j > >m p . The large cross section difference between 
the protons and (p mesons in the hadronic rescattering phase leads to a violation of the 
hydrodynamic mass ordering at low pp i n the final state. 



4. Summary 

We have studied effects of hadronic dissipation on differential elliptic flow in Au+Au 
collisions at RHIC, using a hybrid model which treats the early QGP phase 



Hadronic dissipative effects on transverse dynamics at RHIC 



4 



macroscopically as a perfect fluid and the late hadronic phase microscopically with a 
hadronic cascade. The well-known mass-splitting of the differential elliptic flow observed 
in hydrodynamic models is seen to be mostly generated during the hadronic rescattering 
phase and to be largely due to a redistribution of the momentum anisotropy built up 
during the QGP stage. This redistribution is caused by the mass-dependent flattening 
of the transverse momentum spectra by additional radial flow generated during the 
hadronic stage. The much more weakly interacting mesons do not participate in this 
additional radial flow and thus are not affected by this redistribution of momentum 
anisotropies: their differential elliptic flow remains almost unaffected by hadronic 
rescattering. The net result of dissipative hadronic rescattering is therefore that the 
differential elliptic flow of protons drops below that of the qb mesons, in violation of the 
hydrodynamic mass-ordering. 
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